Draft version August 27, 2012 

Preprint typeset using I^T^^X style emulatcapj v. 5/2/11 



NUCLEOSYNTHESIS IN CORE-COLLAPSE SUPERNOVA EXPLOSIONS 
TRIGGERED BY A QUARK-HADRON PHASE TRANSITION 

NOBUYA NiSHIMURA'^'^, TOBIAS FlSCHER^"', FrIEDRICH-KARL THIELEMANN^'^, 

Carla Frohlich*, Matthias Hempel\ Roger Kappeli\ Gabriel Marti'nez-Pinedo^"*, 
Thomas Rauscher\ Irina Sagert'', and Christian Winteler^ 

e-mail: nobuya.nishimura@unibas.ch 
^ Department of Physics, University of Basel, CH-4056 Basel, Switzerland 
^ GSI, Helmholtzzentrum fiir Schwerionenforschung GmbH, D-64291 Darmstadt, Germany 
Tcchnisclio Univcrsitat Darmstadt, D-64289 Darmstadt, Germany 
Department of Physics, North Carolina State University, NO 27695, USA 
^ Department of Physics and Astronomy, Michigan State University, MI 48824, USA 
Draft version August 27, 2012 

ABSTRACT 

We explore heavy-element nucleosynthesis in the explosion of massive stars which are triggered by 
a quark-hadron phase transition during the early post-bounce phase of core-collapse supernovae. The 
present study is based on general-relativistic radiation hydrodynamics simulations with three-flavor 
Boltzmann neutrino transport in spherical symmetry, which utilize a quark-hadron hybrid equation 
of state based on the MIT bag model for strange quark matter. The quark-hadron phase transition 
inside the stellar core forms a shock wave propagating toward the surface of the proto-ncutron star. 
This shock wave results in an explosion and ejects neutron-rich matter from the outer accreted layers of 
the proto-neutron star. Later, during the cooling phase, the proto-ncutron star develops a proton-rich 
neutrino-driven wind. We present a detailed analysis of the nucleosynthesis outcome in both neutron- 
rich and proton-rich ejecta and compare our integrated nucleosynthesis with observations of the solar 
system and metal-poor stars. For our standard scenario we find that a "weak" r-proccss occurs and 
elements up to the second peak {A ~ 130) are successfully synthesized. Furthermore, uncertainties 
in the explosion dynamics could barely allow to obtain the strong r-process which produces heavier 
isotopes including the third peak {A ~ 195) and actinide elements. 

Subject headings: dense matter - nuclear reactions, nucleosynthesis, abundances - stars: neutron — 
supernovae: general 



1. INTRODUCTION 

Nucleosynthesis from core-coUapsc supernovae is 
certainly responsible for the production of interme- 
diate mass elements including the so-called alpha 
elements and a certain fraction of iron isotopes and 
their neighbors (the Fe-group nuclei). The lighter 
alpha elements from oxygen to silicon have contri- 
butions from hydrostatic burning which takes place 
during stellar evolutions, while the heavier ones 
and the Fe-group isotopes origin at e from explosive 
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related to the source of heavy elements beyond iron. 

There have been strong expectations that the inner- 
most layer of ejecta, close to the forming neutron star, re- 
mains neutron-rich and is a possible site for the r-process. 
For many years the late neutrino-driven wind (NDW), 
following the actual explosion, seemed an adequate r- 
process site in core-collapse supernovae that formed 
a neutron star (lOian fc Wooslevl[l996t iTakahashi et al.l 
119941: IWooslev et al.lll994[ ). Such investigations were fol- 
lowed up by many parametrized calculations that ex- 
plored the sensitivity of the most relevant r-process pa- 
rameter, the neutron-to-seed ratio, via variations of en- 
tropy S*, electron fraction (or number of proton per 
nuclcon) Kg and expansion timescale t (jHoffman et al.l 



(|Thompson et al.ll2001l : iWanajo 1 120061 ) showed that it is 
very hard to attain the required entropies. These results 
have been r ecently confirmed b y fully hydrodynamical 
simulations (jArcones et al.ll2007f ) that in addition showed 
that the presence of a reverse shock does not have any 
major impact on the neutron-to-seed ratio. 

Furthermore, recent investigations noticed that the 
early NDW turns matter proton-rich, producing spe- 
cific Fe-group isotopes and in the subsequent vp- 

process n uclei with a mas s number up to A 2;' 
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200I) . While it was initially hoped that there ex- 
ists a chance that the late NDW still turns neutron- 
rich (after its initial, early proton-rich phase) exist- 
ing core-collapse calculations seem to indicate that the 
wind beco mes even more proton-rich in the lon g -term 
evolution (iFischer et all l2010t iHlidepohl eFal] 120101: 
iFischer et al.ll2012D. although very recent investigations 
(|Martmez-Pinedo et al.l [2011 [Roberts fc Reddv I [20TI 
related to the spectral evolution of v^. and in medium 
might change this somewhat. 

It was realized that the best chance to eject neutron- 
rich matter, is to utilize material stemming from the ini- 
tial collapse and compression, where electron captures 
turned it neutron-rich early in the explosion, before neu- 
trino interactions have the chance to convert it into 
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proton-rich matter in the NDW. Electron-capture su- 
pernovae, which explode without a long phase of accre- 
tion onto th e proto-neutrqn star , apparently provide such 
conditions (jWanajo et alJ[2bllf ). However, the Yg's ob- 
tained under such conditions do not support a strong r- 
process, which successfully reproduces the platinum peak 
of r-elements around A = 195. 

An alternative explosive scenario as a possible site 
for th e r-process has been suggested bv iJaikumar et afl 
(|2007| ). where after a supernova explosion a quark- 
hadron phase transition inside the neutron star can eject 
very neutron-rich material from the neutron star surface. 
This scenario has been referred to as a quark-nova. In 
the present paper, we follow a different approach and 
investigate core-collapse supernova explosions that are 
triggered by a quark-hadron phase transiti on during the 
early post-bounce phase (for details, see iSagert et al.l 
120091 : iFischer et al.ll20Tll ) and their nucleosynthesis fea- 
tures. Under such conditions, zones which are initially 
neutron-rich, can be (promptly) ejected without experi- 
encing the strong effects of the neutrino flux which comes 
from the central proto-neutron star. 

In the following sections, we discuss the supernova ex- 
plosion model and the nuclear physics inputs utilized for 
the nucleosynthesis calculations in Section 2, the condi- 
tions experienced and the resulting ejecta in Section 3, 
and give a detailed analysis of the ejecta composition in 
Section 4, followed by conclusions, also in comparison to 
alternative r-process sites, in Section 5. 

2. METHODOLOGY AND NUCLEAR PHYSICS INPUT 

We explore nucleosynthesis in core-collapse supernovae 
which are triggered by a deconfinement phase transition 
during the early post-bounce phase . The explosion mod - 
els have been discussed in detail bv IFischer et al.l ()20lH ). 
In this section, we summarize the main features of the 
explosion models and the input physics of the nuclear 
reaction network utilized for nucleosynthesis. 

2.1. Hydrodynamics and Equation of State 

A self-consistent supernova explosion model, which we 
adopt to investigate nucleosynthesis in the present work, 
has been carried out by using the AGILE-BOLTZTRAN 
code. This numerical code is based on general-relativistic 
radiation hydrodynamics in spherical symmetry with an 
adaptive grid an d employs a three-flavor B oltzmann neu- 
trino transport (jLiebendorfer et all 120041 ) with detailed 
microphysics, e.g., a realistic nuclear equation of state 
(EOS), weak processes, and (other) nuclear reactions. 
For the current study, the standa rd weak processes w ere 
considered as listed in Table 1 of IFischer et al.l ()2011[ ). 

The baryonic EOS in supernovae needs in general to 
be known for the following three intrinsically different 
regimes: 

1. At temperatures below 6 GK (~ 0.5 MeV), the 
baryon EOS is dominated by heavy nuclei and their 
abundances, determined by individual nuclear re- 
ac tions which a r e not necessarily in equilibrium. 
In IFischer erall ()2010[ ). we included a nuclear re- 
action network, consisting of 20 nuclei, which per- 
mits to simulate a large domain of the progenitor 
up to the helium-layer and is utilized to predict the 
amount of energy generation by explosive nuclear 
burning. 



2. At temperatures above 6 GK (~ 0.5 MeV), the nu- 
clear statistical equilibrium (N SE) can be applied, 
where the baryonic EOS from iLattimer fc Swestvl 
()1991[ ) and iShen et all (jl998l ) are commonly used 
in supernova simulation studies. 

3. Above the nuclear matter density, the state of mat- 
ter is highly uncertain. There exists the possibility 
of a deconfinement phase transit ion. Theref o re, we 
extended the hadronic EOS by iShen et al.l ()1998l ) 
at high densities and temperatures, making use of 
a quark EOS based on the bag model for strange 
quark matter. For the first-order phase transition 
between hadronic and quark phases, we applied 
Gibbs conditions leading to a mixed phase during 
the transition. This results in a continuous phase 
transformation. Details of the quark-hadron hy- 
brid EOS are discussed in IFischer et al.l ()2011|) . 

Besides the contributions from electrons and positrons 
as well as photons. Coul omb corrections to the E OS are 
added by the method of iTimmes fc Arnettl ()1999[) . 

For the current nucleosynthesis predictions, we select 
the explosion calculations of the 10.8 Mq progenitor 
model, where a quark-hadron hybrid EOS was used with 
an early phase transition to quark matter close to nor- 
mal nuclear i natter density (lab eled E0S2, second line 
of Table 2 in IFischer et alj|2011[ ). We chose this model 
because it has an explosion energy consistent with the 
expected order of magnitude of 10^^ e rg (see the second 
line of Table 3 in IFischer erall 120 lit ). The maximum 
gravitational mass for this EOS is 1.5026 Mq. Though 
it is in agreement with the highest precisely known mass 
of a compact star, the Hulse-Taylor pulsar of 1.44 Mq, 
recent mass limits for the physical EOS are based on the 
mihi second pulsars J19 03-f0327 with M = 1.667 ±0.021 
Mq (jFreire et al.ll2011| ) and J1614-2230 with a hig h mass 
of M = 1.97 ± 0.04 Mq (iDemorest et aLllMO) . The 
inclusion of corrections from the strong interaction cou- 
pling constant can stiffen the qu ark EOS and lead to 
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The physical conditions for a possible quark-hadron 
phase transition in the proto-neu t ron st ar are highly un- 
certain. Therefore, IFischer et al.l (|201lD also constructed 
a quark-hadron hybrid EOS which includes corrections 
from the strong interaction coupling constant. They 
showed that such an EOS, with a maximum mass of 
1.67 M0 and a phase transition to quark matter close 
to n uclear saturat i on de nsity (see E0S3 in table 2 and 
3 in IFischer et al.l 120111 ). leads to a qualitatively simi- 
lar explosion scenario in spherical symmetry as obtained 
in the explored models with several EOSs which have 
lower maximum masses. Note, the different post-bounce 
times for the onset of deconfinement, and hence for the 
onset of the explosion, lead to different proto-neutron 
star structures. In general, a delayed phase transition 
at higher densities translates to a longer mass accretion 
phase (for the same progenitor). It results in a more mas- 
sive protoneutron star with a steeper density gradient at 
its surface and a lower Yg at the onset of the explosion 
triggered by the phase transition. An earlier phase tran- 
sition and less steep density gradient would result in a 
slower expansion and a slightly higher Yg. 
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This has important consequences for the nucleosynthe- 
sis of the ejected material since expansion timescale and 
enclosed mass, and also neutrino luminosities and spec- 
tra, de£end_on_theseaspccts (for a detailed discussion, 
see IFischer et al.ll201lD and the behavior of EOSl and 
E 0S2 in their Figure 17. 

W cissenborn et all (|2011[ ) recently showed that it is 
possible to obtain a quark-hadron hybrid EOS which 
allows for both, a maximum mass larger than 2 Mq and 
a low critical density for the appearance of quark matter. 
Whether such quark-hadron hybrid EOS will result in a 
similar dynamical evolution as discussed in IFischer et al.l 
(|201lD will be examined in future simulations. 

2.2. Explosion Scenario 

The supernova post-bounce evolution is characterized 
by mass accretion causing a continuous rise of the cen- 
tral density. Once the central density exceeds the critical 
density for the onset of deconfinement, the quark-hadron 
phase transition takes place, leading to the appearance 
of a quark-hadron mixed phase. Thereby, quark matter 
appears in the supernova core where the highest den- 
sities are experienced. The timescale for the appear- 
ance of quark matter is given by the timescale for the 
central density to rise. This depends on the progenitor 
model, which determines mass accretion rates, and the 
hadronic EOS. Note that the EOS in the mixed phase 
is significantly softer than in hadronic and pure quark 
phases. It is a consequence of the assumed first-order 
phase transition. Mass accretion from the outer layers of 
the progenitor onto the central supernova core leads to a 
continuous rise of the central enclosed mass. When the 
critical mass (given by the hybrid EOS) of the configura- 
tion is obtained, it becomes gravitationally unstable and 
the supernova core begins to contract. The contraction 
proceeds into a collapse which rises the density and con- 
verts the hadronic core into quark matter at around the 
center. A massive pure quark core forms at the center, 
where the EOS stiffens and the collapse halts, and an 
accretion shock forms. The shock wave propagates out 
of the high-density supernova core, remaining an accre- 
tion front with no matter outflow. Once it reaches the 
outer layers of the central core, where the density drops 
over several orders of magnitude, the accretion front ac- 
celerates and turns into a dynamic shock with matter 
outflow. This moment determines the onset of an explo- 
sion, also for supernova models which would otherwise 
not explode in spherical symmetry, based on explosion 
mechanisms discussed so far. Finally, at distances on 
the order of 100 km the expanding shock wave merges 
with the standing shock from the initial core bounce at 
nuclear densities, which was unaffected by the dynamics 
occurring in the supernova core. 

When the (second) shock reaches the neutrino spheres, 
an additional millisecond neutrino burst is released. It 
appears in all flavors, however dominated by and 
i'^/t'> in contrast to the z^e-deleptonization burst related 
to the early bounce shock propagation across the neu- 
trino spheres between 200 and 500 ms after core bounce. 
This second neutrino burst is of particular interest for 
water-Gherenkov neutrin o detectors, which are m ore sen- 
sitive to Ve than to v^. iDasgupta et al.l (l2010f ) demon- 
strated that the currently operating generation of water- 
Gherenkov neutrino detectors (e.g., Super-Kamiokande 



TABLE 1 

Summary of mass zones and their properties 



Zone No. 


M# (IO-^Mq) 


AA/# {Mq) 


i'e.NSE 


*cj (s) 


001-014 


0.000-0.208 


1.496 X 10^* 


0.20 




015-019 


0.210-0.216 


1.474 X 10-5 


~ 0.55 


1.5-4.0 


020-050 


0.217-0.232 


1.063 X 10-5 


~ 0.33 


~ 0.5 


051-120 


0.250-1.482 


1.786 X 10-" 


0.33 ~ 0.50 


~ 0.5 



NSE (below T = 9 GK). t^y. ejection time after the bounce. 

and IceGubc detectors) can resolve such mi Uisecond neu- 
trino burst of the explosion models by IFischer et al.l 
(|20T1 . 

The matter considered for nucleosynthesis studies of 
heavy elements (see Section 3) belongs originally to the 
inner parts of t he silicon and su l fur la yers of the 10.8 Mq 
progenitor from lWooslev et all (|2002[ ). with temperature 
and electron fraction of 3 GK and ~ 0.5, respectively, 
at 800 - 1000 km from the center on the pre-coUapse 
phase. During the collapse and the explosion, material 
contracts and is heated by the shock to temperatures ex- 
ceeding 100 GK and hence completely dissociates into 
free nucleons. At high densities weak processes, mainly 
electron captures, establish a very low proton-to-baryon 
ratio Ye — 0.1 during the core-collapse. The further evo- 
lution is discussed in detail in Section 3. 

2.3. Nuclear Reaction Network 

The nuclear reaction network utilized for the following 
nucleosynthesis simulations of the ejecta is an extension 
of previou s ones, w hich have already been described in 
detail (see iNishimura et al. 2006; Fujimoto et al. 2008). 
The network includes more than 4000 nuclei from neu- 
trons and protons up to fermiu m with atomic number 
Z = 100 (for detail, see Table 1 in INishimura et al.|[2006l) 
and includes proton-rich isotopes as well as neutron- 
rich ones far from stability. It includes two- and 
three-body reactions, decay channels, and electron as 
well as positron cap ture (for details, see network A in 
iFuiimoto et al.|[2007[ ) and screening effects for all rele- 
vant charged- particle reactions. Exp erimentally deter- 
mined masses (|Audi fc Wapstralll995l ) and reaction rates 
are adopted if available. Otherwise, theoretical predic- 
tions for nuclear masses, reaction rates, and beta-decays 
are applied, based o n the finite range droplet mass model 
(|Moller et al.lll995D. Spontaneous and beta-delayed fis- 



sion processes (fStaudt fc Klapdor-KleingrothausI Il992f) 
are taken into account. We adopt the empirica l forrn ula 
for fission fragments bv lKodama fc Ta kahashii (jl975D . 

We also employ neutrino interactions with matter in 
order to include dominant weak interactions affecting 
the evolution of the overall proton/nucleon ratio Y^. For 
(anti-)electron neutrino ca ptures by nucleon s , we a dopt 
reaction rates derived bv iQian fc WooslevI (|1996[ ) but 
ignore any reactions with heavy isotopes because the 
amount of neutrino capture is negligible for the early 
phase of nucleosynthesis. The neutrino fluxes, result- 
ing from the detailed neutrino-radiation hydrodynam- 
ics calculation described in the previous subsection, are 
utilized to determine the actual rates as a function of 
time. These reaction rates depend on the distance from 



4 



Nishimura et al. 




mass zone, # 

Fig. 1. — Initial distribution of mass and Ye as a function of mass 
zone number. Thick lines and dashed lines relate to Lagrangian 
mass coordinate and electron fraction, respectively. Masses are 
measured from the surface of the proto-neutron star (starting at 
for zone 001 of Table [TJ and the electron fractions are adopted at 
the time when the temperature decreases down to T = 9 X 10^ K 
for ejected matter and the end of the hydrodynamic simulation for 
inner non-ejected zones). The plot shows only the mass zone range 
001 to 090 for enclosed mass and Ye. For the mass zones 90 to 
120, Ye stays essentially constant at ~ 0.5 and the enclosed mass 
continues to be proportional to mass zone numbers. 

the proto-neutron star and the mean energy and lumi- 
nosity of neutrinos emitted from the proto-neutron star. 
It depends on the structure and the evolution, which is 
sensitive to the EOS, and the precise evolution history 
of ejected matter including the early phase of the core 
bounce. Thus, this is different from the treatment of 
other nuclear reaction rates, which are determined only 
by local thermodynamic conditions and density and tem- 
perature. 

3. NUCLEOSYNTHESIS IN THE EJECTA 

3.1. Dynamic Evolution of Mass Zones 

In order to calculate the nucleosynthesis evolution 
of ejected matter within a postprocessing approach, 
the dynamic evolution is required in radial Lagrangian 
mass zones. For this reason, the evolution deter- 
mined with the radi ation hydrodynamics code AGI LE- 
BOLTZTRAN (see iMezzacappa &: BruennI 
iLiebendorfer et all 1200 laO bf) which is based on an adap- 
tive grid, was mapped on a Lagrangian grid of 120 mass 
zones. This provides the Lagrangian evolution of phys- 
ical quantities, such as density, temperature, electron 
fraction, and velocity of the ejected material, and in ad- 
dition the neutrino fluxes experienced as a function of 
time. 

The mass zones which are ejected in the explosion arc 
classified in three different categories, related to their 
ejection process and thermodynamic quantities. As listed 
in Table [U zones 001 to 120, given with the ejection 
timescale after bounce {te\) and the final Ye.NSE which 
is Ye at the end of NSE (below T = 9 GK), cover the 
material from the surface of the inner core at 1.48000 
Mq to layers with a corresponding mass of 1.49482 Mq. 
They are also shown with respect to their mass as well 
as Ye-distribution in Figure [T] 

These zones coincide with the matter discussed at the 
end of Section 2.2, where at high densities a Ye decrease 
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Fig. 2. — Top: radial trajectories of mass elements as a function of 
time after bounce. The colors indicate the properties of these mass 
elements: black, red, green and blue lines refer to matter which is 
either (black) not ejected, (red) part of the neutrino-driven wind, 
(green) initially stalled matter which gets boosted by the wind 
and (blue) matter which experiences a prompt ejection. Bottom: 
evolution of Ye as a function of time after the core bounce. The de- 
confinement phase transition happens about 0.4 s after this initial 
bounce, causing the explosion and ejection of matter. The colors 
are the same as in the top panel. 



down to even 0.02 has been noticed. The evolution of 
these mass zones in time is displayed in Figure [2] Zones 
001 to 014 are not ejected within 0.5 s after the core 
bounce. These zones preserve the original low-yg ob- 
tained during collapse and shock wave propagation. As 
they are not ejected, we ignore them in the further nu- 
cleosynthesis discussion, plus all matter originating from 
regions at smaller radii. In Figure [21 they are displayed in 
black. Zones 015 to 019 are ejected in the so-called NDW, 
shown in red. Their is strongly affected by neutrino 
interactions, turning this matter proton-rich. Zones 020 
to 050 (displayed in green) have stalled from infall after 
shock formation and are ejected t hereafter due to neu - 
trino heating and dynamic effects (jFischer et al.l[20Tll ). 
The adjacent zones 051 to 120 are ejected in a prompt 
way, due to the shock wave originating from the decon- 
finement phase transition (displayed in blue). We clearly 
see the division of matter which is ejected in a prompt 
fashion (blue), matter which is coasting and falling in 
again, but gets reaccelerated outward by neutrino en- 
ergy deposition (green) , matter which falls back onto the 
neutron star, but becomes part of the NDW ejecta (red), 
and finally matter which stays on the neutron star and 
will never be ejected (black). 
As shown in the bottom part of Figure [2l the blue and 
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Fig. 3. — Evolution of temperature (top), density (middle), and 
entropy (bottom) of mass zones as a function of time after bounce. 



green zones are neutronized during the collapse via elec- 
tron capture to various degrees, depending on the max- 
imum density attained. Thus the main feature is that 
a strong compression during infall leads to high densi- 
ties (at still low entropies, i.e. highly degenerate matter) 
with large electron Fermi energies which endorse electron 
captures and a strong neutronization of matter with a 
small Ye- Their Ye-values range from 0.35 to 0.5. The 
prompt or quasi-prompt ejection does not change this 
value (with minor effects on the innermost zones, being 
partially affected by the NDW). The effect of neutrino 
and anti-neutrino exposure from the core leads also to an 
increase of Ye, as discussed before. Therefore, the major 
point in favor of an r-process is a fast expansion of that 
material having two aspects: (1) timescales are shorter 
than needed for attaining weak equilibrium and (2) mat- 
ter moves fast to larger radii where the effect of neutri- 
nos vanishes (l/r^). The mass zones in red experience 
a similar effect in their early evolution during collapse, 
but the later evolution leads to values of Yg exceeding 
0.5. This is similar t o recent studies of the NDW (see 
the Introduction and lFischer et alll2010D . where similar 
neutrino and anti-neutrino spectra and fltrx intensities fa- 
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Fig. 4. — Radius of neutrino spheres for electron neutrinos and 
electron anti- neutrinos (top). Luminosities (middle) and mean en- 
ergies (bottom) experienced by mass zone 018, which belongs to 
the neutrino-driven wind with proton-rich ejecta, as a function of 
time after bounce. The difference between the mean energy of anti- 
neutrinos and those of neutrinos is of the order 3 MeV, i.e., less 
than 4A, where A is the neutron-proton mass difference. 



vor proton-rich matter due to the neutron-proton mass 
difference, resulting in different energies available for the 
neutrino/anti- neutrino captures. As the neutrino lumi- 
nosity is still high in the ejection phase, we expect i/p- 
proccss nucleosynthesis. Finally non-ejected mass zones 
(black) can initially also experience interaction with the 
neutrino flux and turn proton-rich while still at larger 
radii and small(cr) densities. Once they settle on the 
surface of the neutron star at high densities, capture of 
degenerate electrons dominates over the neutrino effects, 
and they turn neutron-rich again. 

The thermodynamic conditions (density, temperature, 
and entropy), which are responsible for the nucleosyn- 
thesis results, are shown in Figure [3] for the ejected mass 
zones. These figures indicate a temperature, density, and 
entropy maximum when the quark-hadron phase transi- 
tion occurs, which causes a second core bounce (about 
0.4 s after the first bounce at nuclear densities; see Fig- 



6 



Nishimura et al. 



mass, M# [1 0'^ Mg] 
0.209 0.217 1.482 




mass zone # 

Fig. 5. — Ye and entropy 5, which set the conditions for explo- 
sive nucleosynthesis at the time of matter ejection. The innermost 
ejected zones are proton-rich due the effect of the neutrino-driven 
wind, which also heats matter efficiently, leading to high entropies. 
The outer mass zones, ejected in a more prompt fashion keep 
their original (slightly neutron-rich) from the infall / compression 
phase. 

ure[2]) and an outgoing shock front forms. The expansion 
follows a close to constant entropy, i.e. is adiabatic, once 
matter is ejected (after 0.4 s for the prompt ejection and 
after 1.5 s for the delayed ejection). The matter which 
initially fell back onto the proto-neutron star and is fi- 
nally ejected by the NDW, experiences heating and an 
entropy rise due to this energy deposition by neutrinos. 

3.2. Neutrinos from the Proto-neutron Star 

The properties of the neutrino and anti-neutrino flux 
(luminosities, average energies and neutrino sphere radii) 
can be found in FigureSl It is clearly seen that the initial 
bounce (0 s) at nuclear densities leads to a neutrino burst 
due to electron captures, while the second shock wave 
caused by the quark-hadron phase transition (at 0.4 s) 
also produces antineutrinos. From that point on in time 
the neutrino and anti-neutrino luminosities are compa- 
rable (slightly smaller for anti- neutrinos) . The average 
energies are larger for anti-neutrinos than for neutrinos, 
but the difference remains less than 4 MeV. Neutrino and 
anti- neutrino captures determine the neutron/proton ra- 
tio due to the reactions 

De + p ^ n + i^e + n ^ p + . 

Based on the neutron/proton m ass difference of A = 
1.293 MeV, (jFrohlich et al.ll2006aD could show (see their 
Equation (4)), that with the use of Equations (64a) and 
(64b) in lQian fc WooslevI (fT996h . Ye > in the case that 
the difference between the mean antineutrino and neu- 
trino energies fulfills ep - ej, < 4 A, with e = {E'^)/{E) 
and where (E) is the mean energy and {E'^) is the square 
value of the root-mean-square (rms) energy. Therefore 
Ye > 0.5 is obtained if the timescale for neutrino/anti- 
neutrino captures is shorter than the dynamic timescale. 
Thus, for all conditions discussed here, where neutrino 
and anti- neutrino captures are responsible for the n/p 
ratio, proton-rich conditions are attained, i.e., Yg > 0.5. 
That is exactly what is seen in Figures [2] and |4] for 
mass zones which experience essential neutrino fluxes 
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Fig. 6. — Neutron/seed ratio and remaining mass fractions of 
*He (a-particles) after charged-particle freeze-out. Both properties 
result from the original Ye and entropy S in these mass zones. In 
turn they determine the fraction of heavy elements and whether 
those experience further neutron capture after charged-particle 
freeze-out, which is the key to the pattern of heavy nuclei and the 
maximum mass number attained. In the inner part (green lines) 
an effect is seen, which results from an intermediate fallback before 
final ejection. These mass zones experience first (due to the shock 
from the dcconfinemcnt phase transition) a maximum temperature 
and density, expand afterwards close to adiabatically, heat up dur- 
ing the intermediate fallback, and then expand freely. The line 
indicated with "no-boost" gives the initial neutron/seed ratio after 
the first expansion. The alpha-fraction in the inner zones results 
from the reheating phase. The initial expansion at low Ye's would 
result in vanishing alpha-fractions. 



(weighted by 1/r^) at radii of about 100 km. It is re- 
lated to the spectral evolution of Ve and De after the 
onset of expl o sion and has been discussed in detail in 
iFischer et al.l ()2012l ). For matter at larger radii (about 
1000 km), the timescale for this process is too long and 
minor Ye changes occur, i.e., the initial Ye from the col- 
lapse phase is retained. Matter at smaller radii, on top of 
the neutron star, experiences high densities and electron 
Fermi energies, where electron captures dominate which 
make matter neutron-rich. 

Figure [5] underlines this effect due to neutrino inter- 
actions or electron capture. All outer mass zones keep 
their original Ye, which is due to electron capture at high 
densities during the collapse, and ranges from about 0.48 
(further out) to 0.32 (for the inner quasi-prompt ejected 
matter). Material which fell in initially onto the surface 
of the neutron star and is then ejected via the NDW, has 
been turned proton-rich by neutrino and anti-neutrino 
captures with values up to Ye = 0.55. The NDW also 
leads to energy deposition and an entropy increase to 
maximum values of about 85 ks per baryon. The en- 
tropy in the outer ejected regions is of the order 30-50 
ks per baryon, caused by shock heating during the pas- 
sage of the ejection shock wave. 

3.3. Nucleosynthesis Results 

In the following we show final nucleosynthesis results 
for a number of typical mass zones. In order to get a 
rough idea about the results of explosive nucleosynthesis, 
one can utilize either maximum densities and tempera- 
tures prior to an adiabatic expansion or the entropies at- 
tained in the expanding matter (in radiation-dominated 



NUCLEOSYNTHESIS IN SUPERNOVA EXPLOSIONS TRIGGERED BY A QUARK-HADRON PHASE TRANSITION? 



regimes 5* oc T^/p). 

Co mparing entri es in Figure 5 of ( [Thielemann et aLl 
|1990[) and Fig.3 of ([Thielemann et aLlfl996[ ) leads to the 
conclusion that (1) these are typical conditions for an 
alpha-rich freeze-out from explosive Si-burning and (2) 
one would expect remaining alpha mass-fractions after 
charged-particle freeze-out of the order 20% - 100%. 
One should consider, however, that those calculations 
were performed for hydrodynamic (i.e., free fall) expan- 
sion timescales, which can differ from the actual simu- 
lation, and a value of Ye = 0.4988, i.e. matter neither 
neutron nor proton-rich. Therefore we expect the fol- 
lowing changes: (1) higher /lower entropies within the 
given variety will lead to higher/lower remaining alpha- 
fractions; (2) higher Ke, i-e., more proton-rich matter 
causes an alpha-rich charged-particle freeze-out with re- 
maining free protons, which can thereafter lead to a vp- 
process, if a sufficient flux of electron anti-neutrinos is 
still present; (3) smaller Ye's, i.e., more neutron-rich mat- 
ter permits to bypass the slower triple-alpha reaction via 
the faster aan-reaction, in order to produce heavier nu- 
clei and a reduction in the remaining alpha-fraction is 
expected. In addition, free neutrons are remaining after 
the charged-particle freeze-out. With respect to this lat- 
ter aspect, we expect also the additional behavior: higher 
entropies and lower Fg's lead to a larger neutron/seed 
ratio, seed nuclei being the heaviest nuclei formed after 
charged-particle freeze-out, and permit therefore more 
neutron captures on these seed-nuclei. Dependent on the 
neutron to seed ratio, this could lead to light, medium 
or strong r- processing, producing nuclei in the first, sec- 
ond or third r-process peaks, around A=80, 130 or 195, 
depending on the actual n/sccd ratio attained. 

In Figure [HI the properties of mass zones 020 to 120 
are summarized, those with aYf. < 0.5 which experience 
neutron-rich conditions to a varying degree. Properties 
(1), i.e., the degree of alpha-rich freeze-out and (3), the 
resulting neutron/seed ratio, are displayed. We see a 
complex dependence of these properties on entropy S, Ye, 
expansion timescale r, plus further complications from 
reheating and re-expansion which do not fit to a simple 
expansion interpretation. For similar Ye's, the remain- 
ing alpha-fraction increases with entropy, as expected, 
when looking at the behavior of mass zones 080 to 120. 
Then, following mass zones further in, the decrease in Y,, 
dominates, which permits to pass the alpha-to-carbon 
bottle-neck more efficiently, via the aan-reaction, and 
the remaining alpha-fraction vanishes. As known from 
moderately neutron-rich NDW simulations, none of the 
entropies encountered here leads to sizable neutron/seed 
ratios. For such low entropies, only a strongly neutron- 
rich initial composition permits large (r) neutron/seed ra- 
tios. This is what can be noticed when following mass 
zones from 070 down to 040, where the lowest Ye's are 
encountered. Mass zones 020 to 040 experience a more 
complicated history, initial expansion after shock pas- 
sage, later partial fallback, and then final ejection and 
further expansion. Here, the neutron/seed ratio after the 
first charged-particle freezeout is the one of importance 
for the production of heavy elements. 

The reheating leads to partial photodisintegration 
of heavy elements, the production of alphas, and the 
buildup toward nuclei with mass numbers around A = 



90-110 during the final expansion. The neutron/seed ra- 
tio at this second charged-particle freezeout (boost) is 
rather a measure for local rearrangements of matter. 

The maximum neutron/seed ratio of slightly above 
10 obtained over the complete range of ejected mass 
zones, does not support conditions to produce the third 
r-process peak, in fact only a small production of the 
second peak is expected. The effect of the "boost", i.e. 
second reheating and expansion in mass zones 020-040, 
rearranges/reshapes the abundance distribution via pho- 
todisintegrations and captures, but does not alter this 
conclusion. 

With this background we have a look at the final abun- 
dances of representative mass zones, displayed in Figure 
[7l We sec in the outer layers remaining alpha-fractions 
beyond 60%. These are regions, which experience en- 
tropies of S* = 33-50fcf,/baryon and Ye's close to 0.5 and 
see vanishing neutron/seed ratios. Thus, we expect es- 
sentially the production of the Fe-group up to A = 50- 
70. This can be observed in the last two subfigures of 
Figure [7] (zones 080 and 120). Mass zones 60 and 70, 
which experience the highest entropies and moderately 
decreased Yg's, can move matter up to and (slightly) be- 
yond A — (for mass zone 70 still with a large frac- 
tion of matter in the Fe-group), mostly due to a more 
neutron-rich (lower Ye) charged-particle freeze-out and 
not due to further neutron processing (see neutron/seed 
ratio in Figure IB]). One can also see a variation in the 
final carbon-fraction, underlining how effective the bridg- 
ing of the alpha-to-carbon bottle neck of reactions is in 
comparison to subsequent capture reactions to heavier 
nuclei. At smaller radii (mass zones 020 to 051 and Ye's 
as small as 0.33), also the A = 130 peak starts to be pop- 
ulated. This is due to the neutron/seed ratio of up to 10 
attained after charged-particle freezeout. For the mass 
zones 20-45, abundances are shown after the first expan- 
sion (dashed) and after reheating/boost (solid), which 
is characterized by some photo-disintegration of heavy 
nuclei, the appearance of ^Hc, and further processing of 
nuclei beyond the Fe-group (see also Figure |S]) . Summa- 
rizing the results of all mass zones with Ye < 0.5, we note 
that none of these zones produce matter beyond A ^ 130, 
nor show the 130 peak dominating abundances. 

Deeper mass zones, ejected via the NDW turns proton- 
rich and they experience entropies as high as 85fcs 
per baryon. These are conditions where we expect 
a i/jp-process (jFrohlich et all I2006bl: IPruet at al.l [200l 
IWanaiol[2006|) . This can be seen in the first three pan- 
els of Figure [71 dashed lines show abundances without 
the inclusion of the i/p-process, solid lines show the fi- 
nal results after fp-processing. While in terms of total 
mass fraction the production of nuclei beyond the Fe- 
group is not too impressive. Figure [3] displays this more 
prominently, where the overproduction ratio over solar 
is plotted. In the proton-rich environment anti-neutrino 
captures on free protons produce neutrons and permit to 
overcome the rp-process waiting point ^"^Ge via an (n,p)- 
rcaction, winning against a slower /3+-decay. In this way, 
nuclei up to A = 80 — 90 can be produced on the proton- 
rich side of stability. The comparison of the black and 
red bullets indicates the strength of this process. 

4. SURVEY AND INTEGRATED COMPOSITION 
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Fig. 7. — Final nucleosynthesis results for selected mass zones showing mass fractions as a function of nuclear mass number A. The mass 
zone number is given in each panel. In panels, where wc show dashed and solid lines, the solid lines correspond always to the final result. In 
zones up to 19 the dashed lines correspond to abundances without including the !^p-process, i.e. before proceeding beyond the beta-decay 
bottle-neck "^"Ge. In zones 20-45 they show abundances before the reheating boost. 



After having discussed the individual composition, 
ejected from different positions in the exploding model, 
we want to give a final survey of the conditions attained 
in all mass zones of explosive Si-burning which are af- 
fected by the explosion mechanism (here the deconfine- 
ment phase transition). The discussion of nucleosynthe- 
sis in layers further out is only affected by the energy in 
the shock wave, as has been discus sed extensively in the 
literature ([Thielemann et aTlllQQGf) . and will not be re- 
peated here. Finally, we also discuss the overall features 
of the integrated yields (of these inner mass zones, close 
to the explosion mechanism). The survey is displayed 
in Figure [TU] which features abundances after charged- 
particle freeze-out, and thus the setting for the final nu- 
cleosynthesis features, if free neutrons or protons are re- 
maining in sizable fractions. 

Essentially, all mass zones have total entropies in ex- 
cess of S' = SOfcs per baryon. In the outer mass zones 
with a Ye close to 0.5 this produces high alpha-fractions 
plus dominantly ^^Ni. Moving somewhat further in with 
slightly increasing entropies, noticeable amounts of ^"^Ge 
are produced as well, with the ^''^Ge/^^Ni being a measure 
of entropy. Decreasing Yg, i.e., having more neutron- rich 
conditions, changes the initially pure alpha-rich charged- 
particle freezeout to an alpha-rich freeze-out with siz- 



able abundances up to A = 90. A further decrease in 
Yg causes remaining free neutrons, which permit an ad- 
ditional sequence of neutron captures. The decrease of 
Ye down to 0.33, leads to a charged-particle freeze-out 
with vanishing alpha-fractions but a sufficient amount of 
free neutrons, which permit later to produce nucleosyn- 
thesis ejecta in the mass A=130 peak. The mass zones 
which are affected by a reheating boost are characterized 
by photo-disintegrations and a second charged-particle 
freeze-out with remaining alpha-fractions, Finally the in- 
nermost proton-rich zones, with higher entropies of about 
S = Sbks per baryon and experiencing a continuous 
neutrino-fiux, show an alpha-rich and proton-rich freeze- 
out with Yg's up to 0.55. This causes a fp-process after 
charged-particle freeze-out , permitted by neutron pro- 
duction via anti-neutrino captures on free protons. This 
process will produce nuclei up to A = 90, but on the 
proton-rich side of stability, especially Sr, Y, and Zr iso- 
topes. 

After this survey of the conditions at charged-particle 
freeze-out, combined with the final ejecta composition 
as a function of radial mass coordinate as discussed in 
the previous section, we want to give an integrated pre- 
sentation for these inner mass zones which experience 
conditions for the possible formation of nuclei beyond 
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Fig. 8. — Top: The comparison of results with (solid) and with- 
out (dashed) the effect of reheating and second expansion (boost) 
integrated over the whole range of mass zones 20-50. Bottom: The 
ratios of the two cases, based on the full hydrodynamic evolution 
with reheating (boost) and the neglcction of reheating around 1.5 
s.The reheating leads to photodisintcgrations, a strong appearance 
of *He and a reshaping of the heavy element distribution. 
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Fig. 9. — Resulting nucleosynthesis for the entire mass zones in 
the neutrino-driven wind, experiencing a i/p-process. The result 
is shown for two options; (a) including the anti-neutrino captures 
which turn protons into neutrons in the late phase of nucleosynthe- 
sis and permit to overcome the ^*Ge beta-decay bottle neck via an 
(n,p)-reaction (filled red circles), (b) neglecting this effect (open 
black circles). It can be seen that the inclusion of anti- neutrino 
reactions enhances abundances of nuclei heavier than A = 64 and 
permits the production of nuclei up to A = 80 — 90. 



the Fe-group. Figure [TT] (top) shows the composition 
for this range in mass numbers in comparison to solar 
r-abundances. As we do not yet know the frequency of 
such events, i.e. which range of the initial mass function 
of stellar masses leads to these types of explosions, we 
show a scaling normalized to the A ^ 100 mass region 
(where the dominant abundances are obtained). The 
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Fig. 10. — Overview of dominating conditions, i.e. abundances 
of a few key nuclei after chargod-particle freeze-out, as a func- 
tion of radial Lagrangian mass, with variations from (a) typical 
explosive Si-burning products and an alpha-rich freeze-out in the 
outer zones to (b) an increasing remaining neutron abundance after 
charged-particle freeze-out, permitting a weak r-proccss, down to 
(c) proton-rich neutrino-driven wind ejecta, permitting the onset 
of an i/p- process. 



bottom part of Figure [TT] shows the individual contri- 
butions to the overall abundances by the different mass 
zones as presented in Table [T] and Figure [TJ and shown 
in different colors in Figure [2] and [3] as well as [SJ |6l and 
[71 The high end of the abundance distribution in the 
prompt (blue, dashed) as well as the delayed (boosted) 
ejecta (green, solid) are due to the most neutron-rich 
conditions with close to 0.33. The outermost ejected 
mass zones with Ye closer to 0.5 contribute to the Fe- 
group and matter up to A = 80 (lower range of mass 
numbers of blue dashed line) . The NDW (red line) with 
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Fig. 11. — Integrated abundance distributions as a function of 
atomic mass number A for all ejccta in comparison to solar r- 
abundances (top) (normalized to the A=100 region) and separated 
by ejection process ('feoiiomj. We also indicate the effect of uncer- 
tainties in the Ye determination, given in terms of percentage Pcor, 
discussed in section 5. 
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Fig. 12. — Density evolution as a function of time. Shown is the 
standard case with the original explosion simulation according to 
EOS2 and exponential density expansions fitted to this case (/ = 1) 
as well as slower expansions with longer expansion timescalos by a 
factor 2, 5, and 10. 



slightly proton-rich ejecta { — 0.55) produces signifi- 
cant abundances only up to ^ = 64, if the same normal- 
ization for the relevant mass zones is used (see, however, 
also Figure O ■ 
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Fig. 13. — The corresponding temperature evolutions with adia- 
batic index 4/3. 



It is clearly visible that, first of all, objects like these su- 
pernovae, exploding by a mechanism based on EOSs with 
a low-density quark-hadron phase transition, do only ex- 
perience a weak r-process in ejected mass zones which 
were neutronized during collapse. There is no matter 
produced in the third r-process peak. If normalizing the 
abundance curve at A = 100, in order to avoid an over- 
production of this mass region, also only a small contri- 
bution (less than 10% is expected to the second r-process 
peak, i.e., A = 130). The major production affects the 
atomic mass range from A = 80 to 115, curiously also 
reproducing a minimum at ^ = 97-99. Thus, the type 
of events discussed here, contribute to the whole mass re- 
gion beyond the Fe-group up to A = 115 in a significant 
way, accompanied by minor contributions to the second 
r-process peak a.t A = 130. 



5. MODEL UNCERTAINTIES 



5.1. 



Equation of State Uncertainties 

The major uncertainty in the present calculations is 
expected from the uncertain properties of the EOS, es- 
pecially the quark-hadron phase transition which is at 
the origin of the explosion mechanism discussed here. 
This was already revi ewed extensiv e ly in Section 2.1. 
The analysis done by iFischer et al.l ()201lD with varia- 
tions in the properties of the phase transition, led to dif- 
ferent post-bounce times for the onset of the explosion 
and different neutron star structures. Later transitions 
at higher densities cause more massive progenitors with 
steeper density gradients at the surface. The standard 
case of th e present pap e r corr esponds to E0S2 in Fig- 
ure 17 of IFischer et all (poTTI ). i.e., a late(r) explosion 
with steep(er) density gradient and thus a fast expansion. 
Variations toward the properties of EOS 1, with an earlier 
explosion, would cause slower expansions. Therefore, we 
performed test (nucleosynthesis) calculations with slower 
expansions for a specific mass zone (here zone 51, see 
also Figures [7]) which provided among the best r-process 
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conditions in our standard case. FigslT^ and [T3] show 
variations with exponential density expansions fitted to 
the standard case (/ = 1) and slower expansions with 
longer expansion times scales by a factor of two, five, 
and !0 (choosing the temperature according to adiabatic 
expansions with the same entropy). 
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Fig. 14. — Final abundances for mass zone #51 with differ- 
ent expansion speeds. The slower expansions (/ = 2, 10) lead 
to smaller remaining alpha fractions (as expected), but also to a 
slightly stronger (weak) r-process, indicated by the abundances in 
the A = 130 peak. 

In typical NDW environments for given S and Y^, a 
fast(er) expansion leads to a strong(er) alpha-rich freeze- 
out, thus more alpha particles and less seeds and a higher 
neutron/seed ratio, promoting the options for a stronger 
r-process. However, in the mass zones of interest here, 
the Ke is quite low and the remaining alpha-fraction 
close to negligible. This feature was discussed exten- 
sively in Section 3.3 with respect to Figures |6] and [T] 
We see the effect in Figure [Ml that the fastest expan- 
sion {f — 1) shows the highest remaining alpha- fraction, 
but only of the order 10~^. This behavior lowers the 
initial seed abundance after chargcd-particle freeze-out 
by the same amount, which is negligible. This causes 
essentially the same seed abundances and neutron/seed 
ratio for all calculations with different expansion speeds 
(/ = ! - 10) and cancels this dominant effect in typ- 
ical neutrino-driven wind environments. In fact, the 
slowest expansion leads to the highest abundances at 
A = !30. Here, a second-order effect seems to take 
over, the different density and temperature conditions 
for slower expansions lead to an apparently slightly more 
neutron-rich r-process path, encountering smaller beta- 
decay half-lives which permit to proceed faster to heavier 
nuclei. This behavior, as shown in Figure I14i indicates 
a slightly stronger (weak) r-process for the slower ex- 
pansions. However, the effect is not strong, only clearly 
visible for / = !0, which is larger than the uncertainty 
encountered in our EOS variations. 

5.2. Ye Uncertainties 

One can argue, that there might exist some uncertainty 
in weak interactions (el ectron/positron captures and 
neutr ino captures; sec, e. g.lLanganke fc Martmez-Pinedol 
[2003t IJanka et al, ,2007HCole et al.ll20!2[) which deter- 



mine Kg, especially due to recent developments of in- 
cluding charged-current r ates that are consistent with 
the EOS (for details, s ee iMartinez-Pinedo et al.l I20!2l : 
[Roberts fc Reddv"1l20!2l ). These very recent preliminary 
investigations include medium effects for the (electron) 
neutrino and antineutrino capture reactions, which lead 
to changed reaction Q-values by adding nucleon interac- 
tion potential differences for neutrons and protons, en- 
hancing as a result the difference between average anti- 
neutrino and neutrino energies. This causes a change in 
Yf,. All our earlier discussions had the main emphasis 
that the influence of the neutrino fiux from the proto- 
neutron star would cause an incr ease in Y^. The prelim- 
inary analy sis of such effects by IMartinez-Pinedo et al.l 
(|2012D and [Roberts fc Reddv I ([201^, which were not 
yet included in the present calculations, shows that the 
Ye increase is weakened or even moderately neutron-rich 
conditions (Ye ^ 0.5) can result. 

For this reason, we also repeated the present nucleosyn- 
thesis calculations with variations in the initial Ye for the 
mass zones experiencing prompt and delayed explosions, 
according to the following recipe: 

Ye^ cor = 0.5 + (Ye - 0.5) X (f + ^) , 

where the Ye, cor's are the corrected ones and Pcor denotes 
the percentage of uncertainty in deviations of Ye from the 
symmetric value 0.5, which enlarges these deviation from 
0.5. This lowers the initially only moderately neutron- 
rich Ye's in regions which already produced A = 130 
nuclei. The nucleosynthesis results are presented in Fig- 
ures [TT] and [15] and show the options of obtaining a full or 
weak r-process. However, we expect that any uncertain- 
ties beyond 20% are unrealistic for the explosion model 
we adopt in the current work. 

As is obvious from the discussion above, that only a 
"weak" r-process can be supported by the nucleosyn- 
thesis conditions found in the explosion mechanism dis- 
cussed and presented here, one might wonder whether 
such conditions support abundance features found in 
"weak r-pro c ess" l ow-metallicity stars as observed by 
[Honda et all ([2006[ ). For this reason, we also show such 
a comparison in Figure [T5l What can be seen is that 
these observations also show sizable r-process features 
above the A = 130 peak, although weaker than in so- 
lar r-element abundances. If such abundance distribu- 
tions are the result of a single nucleosynthesis pollution, 
also the "weak" r-process found in the present paper can 
only marginally explain such features. One could argue, 
however, that such observed abundance features are a 
combination of at least two pollutions, one (low level) 
solar r-contribution plus another "weak" r-contribution 
extending only up to A = 130. 

5.3. Multi- Dimensional Effects 

Nucleosynthesis uncertainties due to multi-dimensional 
effects, i.e., convective turnover, have not been consid- 
ered here, but we want to discuss their possible impact. 
Note that convection and the presence of unstable fluid 
motion are typically associated with compositional mix- 
ing, also of the ejected material. Their development, 
e.g., after the shock passage across a low-density envi- 
ronment, can le ave peculiar features such as discussed in 
([Wanaio fc Janka.2012, ) at the example of a S.SM© star 
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explosion. In order to give an estimate for our nucleosyn- 
thesis results obtained in the presence o f possible convec- 
tion, we evaluated the Ledoux criterion (jWilson fc Mavld 
119881 ). It is related to entropy per baryon and lepton- 
number gradients. We find that regions where these gra- 
dients become largely negative, i.e. regions where con- 
vection may develop, correspond to either fall-back of 
material or the shock propagation in Yg ~ 0.5 mate- 
rial. Thus, matter with such uncertainties is either not 
ejected (fallback) or mixing occurs in regions with very 
similar Yg. Hence, we expect the impact of possible con- 
vection on the here discussed nucleosynthesis scenario to 
be of minor importance, in particular the ejected mate- 
rial which attains low-Fg is not affected. Independent 
of this expected unimportance of convection for nucle- 
osynthesis results, we should keep in mind that matter 
is strongly accelerated and thus ejected very fast. Even 
if instabilities exist, they should not lead to significant 
mixing on such small timescales after the onset of the 
explosion and affect the nucleosynthesis. 
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Fig. 15. — Integrated abundance distributions as a function of 
atomic charge number Z in comparison to abundan ce features in 
"weak r-process" low metallicity stars as observed bv lHonda et al.l 
II2006I ). We also indicate the effect of uncertainties in the Ye deter- 
mination, given in terms of percentage Pcor 



Nevertheless, this nucleosynthesis site inhabits another 
possibility of convection due to the explosion mechanism 
related to the quark-hadron phase transition. It takes 
place on timescales on the order of several 100 ms, de- 
pending on details of the quark-hadron EOS, in the ab- 
sence of any other, e.g., neutrino-driven explosion. On 
these timescales, large-scale convection may develop be- 
fore the transition/explosion is initiated. Exploring its 
impact on the entire explosion and nucleosynthesis sce- 
nario will require multi-dimensional simulations, espe- 
cially related to the point in time when this explosion 
starts (possibly affected by convection and related neu- 
trino losses and/or uncertainties in the EOS, which both 
determine the Ye at the onset of the explosion) . This is 
beyond the scope of the present article and will require 
future investigations. 

6. DISCUSSION AND SUMMARY 

Supernova nucleosynthesis is well understood for 
the outer ejected mass zones, which can be well 
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Thielemann et al.l 1201 ID . What remains uncertain is 
the composition of the innermost ejecta, directly linked 
to the explosion mechanism, i.e., the collapse and 
explosion phase. In the present paper, we analyzed 
these mass zones of core-collapse supernova explosions 
triggered by a quark-hadron pha se transition dur- 
ing the early pos t-bounce phase (jSagert et al.l 120091 : 
iFischer et aL|[20Tll ). A number of aspects are important 
for understanding these results. The very innermost 
ejecta are strongly affected by the NDW. Recent 
investigations noticed that this NDW turns matter 
proton-rich, producing specific Fe-group isotopes and 
in the subsequent i/p-process nuclei with masses up 
to A = 80 - 90 jL icbendorfcr et al.l 120031: iPruet at al " 



2n05| :lFrohlich et aT.li2006a.b : Pruet at al.ll2006l: iWanaio 



2006f ) . Even in t he long-term evolution proton- rich 
conditions prevail (IFischer et al.l 120101 : iHiidepohl et al.l 
l20T0f ). Thus, there seems to exist no chance to produce 
r-process matter in these innermost regions, despite 
many interesting parameter studies for NDW ejecta in 
terms of ent ropy S, electron fraction Y^, and expansion 
timcscalc t (Hoffman et al."1997^ 'M ever fc Brownlll997l: 
l^iburghaus ct al. 1999a; Farouqi et al.l I2010() ~or hy- 
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mm and 



explored charged-current 
weak processes (i.e., electron- neutrino and anti- neutrino 
captures on neutrons and protons), consistent with the 
EOSs in neutrino-driven supernova explosion. They find 
modifications of the reaction Q-values due to medium 
effects, which increase spectral differences between 
and v^. This results in slightly neutron-rich ejecta. 
Nevertheless, since the entropy per baryon is still low a 
strong 7'-proccss is unlikely to occur. 

In order to obtain r-process conditions, a better chance 
to eject neutron-rich matter is provided, when neutron- 
rich matter stems from the initial collapse and com- 
pression, where electron captures made it neutron-rich, 
early in an explosion, before neutrino interactions have 
the chance to turn it proton-rich in the NDW. Core- 
collapse supernovae, exploding via the quark-hadron 
phase transition (the focus of the present stu dy), or 
electron-capture supernovae (jWanaio et al.|[20Tl[) , which 
explode without a long phase of accretion onto the proto- 
neutron star, both lead to a rather prompt ejection of 
prior compressed and neutronized matter. However, the 
Yg obtained under such conditions does not support a 
full r-process. This kind of outcome also characterizes 
the conditions we find in the prompt and quasi-prompt 
ejecta of the present study, which did not experience a 
strong neutrino flux. However, the Yg-values attained are 
not smaller than 0.33. Whether uncertainties in the in- 
put or explosion physics can change this down to values 
close to 0.23, necessary for obtaining the third r-process 
peak, remains to be shown in an upcoming article where 
we will include charged-current reaction rates th at are 
consistent with the EOS following , Martinez-Pinedo et al.l 
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((Ml) and [Roberts fc Reddv I (IMl) . 

A rela ted, but different, ph enomenon has been dis- 
cussed in lJaikumar et all (|2007l , and references therein). 
In their scenario, the quark-hadron phase transition 
does not happen shortly after core collapse and is not 
the cause of the supernova explosion. Instead a regu- 
lar neutrino-driven supernova explosion occurs first and 
leads to a deleptonized neutron star with low-Fe and a 
steep density gradient at the surface. Their scenario of a 
"quark-novae" , occurring via a quark-hadron phase tran- 
sition in an existing neutron star, is reported to produce 
r-process elements via the prompt ejection of extremely 
neutron-rich material from the neutron star surface. In 
this sense, quark-novae differ from the scenario explored 
in the current study. They result in a more compact 
object already in /3-equilibrium, such that the expansion 
timescale of the ejected material is much faster, com- 
pared to th e scenario discussed i n the present paper. 
However, in iJaikumar et all ()2007[) the authors apply a 
very simplified description for the quark-hadron phase 
transition and ignore neutrino transport in general, as 
well as weak processes in the hadronic part of their neu- 
tron star. The latter may be essential for producing neu- 
trinos which may change Ye during shock passage and 
hence lead to differen t nucleosynthesis result s than what 
has been discussed in lJaikumar et ahl (120071 ). It remains 
to be shown if the obtained favorable conditions for the 
r-process will remain, applying more sophisticated input 
physics. 

Given our results described above, we conclude that su- 
pernova explosions triggered via the quark-hadron phase 
transition during the early post-bounce phase can con- 
tribute to a weak r-proce ss, consistent with o bservations 
of low-metallicity stars fHonda et al.' '2006|) and with 
LEPP abundances (Travaglio et al. 2004). Such condi- 
tions apparently do not occur in regular core-collapse 
supernovae. However, strong r-process conditions, which 
also produce the third r-process peak and the ac- 
tinides, have in simulations only materialized in neutron 
star mergers (Freiburghausi et all Il999bl : iGorielv et all 
120 lit ) , fast rotating corc-collapse s upernovae with 
strong magnetic fields and jet ejecta (iCanrero ^ [200l 
iNishimura et all l2006t iFuiimoto etldl I2008D or accre- 
tion disks arouiid bl ack holes (jSurman et al.l 120081 : 
iWanaio fc Jankal |2012[ ). For the explosion model dis- 
cussed here, investigation of uncertainties in the hydro- 
dynamics and the nuclear physics input showed that a 
strong r-process to produce up to the third peak ele- 
ments is rarely possible to obtain. 

Neutron star mergers have been shown to be power- 
ful sources of r-process matter, in fact ejecting a fac- 
tor of 100 - 1000 more r-process material than required 
on average from core-collapse supernovae, if those would 



have to explain solar r-process abundances. This would 
actually support the large scatter of Eu/Fe in compari- 
son to, e.g., O, Mg, Si, S, and Ca/Fe, where the latter 
arc clearly produced in supernovae. The only problem 
is that it might be hard to explain the early appear- 
ance of r-process matter for mctallicities at and below 
[Fe/H] = —3. Neutron star mergers will appear after the 
first supernovae hav e already produced Fe and studies by 
lArgast et al.l (|2004D showed that one expects r-process 
products only for metallicities [Fe/H] = —3 ^ —2. Some 
recent studies, which include the fact that our Galaxy 
is possibly the result from smaller merging subsystems 
(with different star formation rates) have been expected 
to show a way out of this dilemma. If this cannot 
solved, we need another strong r-process source already 
at low metallicities, and possibly jets from rotating core 
collapses with st rong magnetic fields could be the so- 
lutioii (see, e.g. |Nishimura et al.l 120061 : iFujimoto et al.l 
l2007tlWinteler et a\\\20lK 
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